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Abstract 

Some difficulties are pointed out in the methods for identification of obstacles based 
on the numerical verification of the inclusion of a function in the range of an operator. 
Numerical examples are given to illustrate theoretical conclusions. Alternative methods 
of identification of obstacles are mentioned: the Support Function Method (SFM) and 
the Modified Rayleigh Conjecture (MRC) method. 
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1 Analysis 

During the last decade there are many papers published, in which methods for identi- 
fication of an obstacle are proposed, which are based on a numerical verification of the 
inclusion of some function / := f{a, z), z G M^, a G 5^, in the range -R(-B) of a certain 
operator B. Examples of such methods include ISIjIH]- It is proved in this paper 
that the methods, proposed in the above papers, have essential difficulties. This also is 
demonstrated by numerical experiments. Although it is true that / ^ -R(-B) when z ^ 
it turns out that in any neighborhood of /, however small, there are elements from R{B). 
Also, although / G R{B) when z & D, there are elements in every neighborhood of /, 



however small, which do not belong to R{B) even if z ^ D. Therefore it is not possible 
to construct a stable numerical method for identification of D based on checking the 
inclusions / ^ R{B) and / G R{B). 

We prove below that the range R{B) is dense in the space L^(S'^). 

Assumption (A): We assume throughout that k'^ is not a Dirichlet eigenvalue of the 
Laplacian in D. 

Let us introduce some notations: N(B) and R(B) are, respectively, the null-space and 
the range of a linear operator S, D G is a bounded domain (obstacle) with a smooth 
boundary S , D' = \ D , uq = e*'^"'^, k = const > 0, a G S"^ is a unit vector, N is the 
unit normal to S pointing into D', g = g{x,y,k) := g{\x — y\) := f^^j^z^, / := e~*'^"''^, 
where ^ G and a' G S^, a' := xr~^, r = \x\, u = u{x, a, k) is the scattering solution: 

(A + k'^)u = m D', u\s = 0, (1) 

u = Uo + v, V = A{a' , a, k)e^'^^r~^ + o{r^^), as r oo, xr~^ = a', (2) 

where A := A{a', a, k) is called the scattering amplitude, corresponding to the obstacle 
D and the Dirichlet boundary condition. Let G = G{x, y, k) be the resolvent kernel of 
the Dirichlet Laplacian in D': 

{A + k^)G = -6{x-y) in D',G\s = 0, (3) 

and G satisfies the outgoing radiation condition. 
If 

{A + k^)w = in D\ w\s = h, (4) 
and w satisfies the radiation condition, then ([7j) one has 

u,W=/'g„(..,)M.M., u, = .(.',*)e«=V-H-o(.-), a. r ^ oo, 

(5) 



s 



We write a(a') for a(a', fc), and 

a{a) := Bh := — / UN{s^—a')h{s)ds^ (6) 
4vr Js 

as follows from Ramm's lemma: 
Lemma 1 (pj, p. 46) One has: 

G{x,y, k) = g(r)u{y, —a' , k) + o{r^^), as r = \x\ ^ oo, xr^^ = a', (7) 

where u is the scattering solution (l)-(2). 
One can write the scattering amplitude as: 

Aia\a,k) = -— I UN(s,-a')e''"''ds. (8) 
4vr Js 



The following claim is proved in jH]: 

Claim: f := e"^^"''^ e R{B) if and only ii z e D. 

Proof of the claim. Our proof is based on the results in [7]. 

a) Let us assume that / = Bh, i.e., / G R{B), and prove that z E D. Define 
p{y) '■= giU:^) ~ 4'{y)j where i/jIu) := JgGN{s,y)h{s)ds. The function p{y) solves the 
Helmholtz equation (4) in the region \y\ > \z\ and p{y) = o(^) as \y\ oo because of 
(7) and of the relation Bh = f. Therefore (see [7], p. 25) p = in the region \y\ > \z\. 
Since ip is bounded in D' and g{y, z) ^ oo as y z, we get a contradiction unless z E D. 
Thus, / G R{B) implies z E D. 

b) Let us prove that z E D implies / G R{B). Define ipi^y) := JgGN{s,y)g{s, z)ds, 
and h := g{s, z). Then, by Green's formula, one has ip{y) = g{y, z). Taking here \y\ —>■ oo, 
ill = a', and using (7), one gets / = Bh, so / G R{B). The claim is proved. □ 

Consider B : L'^{S) L\S^), and A : L'^{S'^) L\S^), where B is defined in (6) 
and Aq := J^a A{a', a)q{a)da. 

Theorem 1. The ranges R{B) and R{A) are dense in L^(S'^). 

Proof. Recall that assumption (A) holds. It is sufficient to prove that N{B*) = {0} 
and N{A*) = {0}. Assume = B*q = J^^ un{s, —a')qda', where the overline stands 
for complex conjugate. Taking complex conjugate and denoting g by g again, one gets 
= Jg2UN{s, —a')qda' . Define w{x) := Jg2u{x, —a')qda'. Then w = w^^ = on 
S, and w solves equation (1) in D'. By the uniqueness of the solution to the Cauchy 
problem, w = in D'. Let us derive from this that q = 0. One has w = wq + V, where 

:= '^Qda' , and V := jg2v{x, —a' ,k)qda' satisfies the radiation condition. 

Therefore, wq{x) = in D\ as follows from Lemma 2 proved below. By the unique 
continuation, wo{x) = in M^, and this implies g = by the injectivity of the Fourier 
transform. This proves the first statement of Theorem 1. Its second statement is proved 
below. □ 

Let us now prove Lemma 2, mentioned above. 
We keep the notations used in the above proof. 
Lemma 2. li w = wq + V = Q in D', then = in D'. 

Proof. The idea of the proof is simple: since wq does not satisfy the radiation 
condition, and V satisfies it, one concludes that wq = 0. Let us give the details. The key 
formula is ([7j, p. 54): 

/ e''='^-%(/?)d/3 = — [7g(-«)-7g(«)] + o(-), r-.oo, (9) 

where 7 := e*'^'^/r, and one assumes q G C^^S"^). 

If r := |x| —>■ 00, then, by Lemma 2, assuming q G (7^(5*^), and using the relation 
w = Wq + V = in D', one gets q{a) = for all a G S*^. Thus, Lemma 2 is proved 
under the additional assumption q G C^{S^). If g G L^(5'^), then one uses a similar 
argument in a weak sense, i.e., with x := rj3, one considers the inner product in L'^{S'^) 
of wo{r(3) and a smooth test function h G C°°(S'^), and applies Lemma 2 to the function 
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J^^e-'^^'-^'hdp. Then, using arbitrariness of h, one concludes that g = as an element 
of L^(S'^). Lemma 2 is proved. □. 

Let us prove the second statement of Theorem 1. 

Assume now that A*q = 0. Taking complex conjugate, and using the reciprocity 
relation: A{a,P) = A{—P, —a), one gets an equation: 

/ A{a,l3)hd(3 = 0, (10) 

where h = g(— /?). Define w{x) := Jg2u{x, P)hdp. Then w = Wq + V, where Wq : = 
Jg2 G^^'^'^hdp, and V := fg2v{x, P)hdp satisfies the radiation condition. Equation (10) 
implies that = o(^) as r — oo. Since function V solves equation (1) and V = o(^), 
one concludes (see jTj, p. 25), that = in D', so that w = Wq in D'. Thus, Wo\s = 
w\s = 0. Since Wq solves equation (1) in D and Wo\s = 0, one gets, using Assumption 
(A), that Wq = in D . This and the unique continuation property imply wq = in R^. 
Consequently, /i = 0, so g = 0, as claimed. Theorem 1 is proved. □ 

Remark 1. In |2j the 2D inverse obstacle scattering problem is considered. It is 
proposed to solve the equation (1.9) in j2j: 

[ A{a,p)gdp = e-'^''-\ (11) 

where A is the scattering amplitude at a fixed k > 0, is the unit circle, and z is a 
point on M^. If Q = z) is found, the boundary 5* of the obstacle is to be found by 
finding those z for which ||^|| := ||^(/3, 2;) 112,2(51) is maximal. Assuming that k"^ is not a 
Dirichlet or Neumann eigenvalue of the Laplacian in that D is a smooth, bounded, 
simply connected domain, the authors state Theorem 2.1 i;2j, p. 386, which says that for 
every e > there exists a function Q G L^(S'^), such that 

lim||6;(A^)|| = oo, (12) 

and ( see [2], p. 386), 

A{a,l3)gdp-e-'^''-'\\<e. (13) 
There are several questions concerning the proposed method. 

First, equation (11), in general, is not solvable. The authors propose to solve it 
approximately, by a regularization method. The regularization method applies for stable 
solution of solvable ill-posed equations (with exact or noisy data). If equation (11) is not 
solvable, it is not clear what numerical " solution" one seeks by a regularization method. 

Secondly, since the kernel of the integral operator in (11) is smooth, one can always 
find, for any z G M^, infinitely many Q with arbitrary large ||^||, such that (13) holds. 
Therefore it is not clear how and why, using (12), one can find S numerically by the 
proposed method. 
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Remark 2. In [2], p. 386, Theorem 2.1, it is claimed that for every e > and every 
7/0 e -D there exists a function Q such that inequahty (13) (which is (2.8) on p. 386 of [2]) 
holds and 1 1^| | ^ oo as ?/o dD. Such a ^ is used in [2] in a "simple method for solving 
inverse scattering problem". However, in fact there exist infinitely many Q such that 
inequality (13) holds and ||^|| —>■ oo, regardless of where yo is. Therefore it is not clear 
how one can use the method proposed in [2] for solving the inverse scattering problem 
with any degree of confidence in the result. 

Remark 3. In P it is mentioned that the methods (called LSM-linear sampling 
methods) proposed in papers [SI, 0, [HI produce numerically results which are inferior to 
these obtained by linearized Born-type inversion. There is no guarantee of any accuracy 
in recovery of the obstacle by LSM. Therefore it is of interest to experiment numerically 
with other inversion methods. In [7], p. 94, (see also [H] , [HI , [I2| ) a method (SFM-support 
function method) is proposed for recovery of strictly convex obstacles from the scattering 
amplitude. This method allows one to recover the support function of the obstacle, and 
the boundary of the obstacle is obtained from this function explicitly. Error estimates of 
this method are obtained for the case when the data are noisy j7], p. 104. The method is 
asymptotically exact for large wavenumbers, but it works numerically even for fca ~ 1, 
as shown in 131. For the Dirichlet, Neumann and Robin boundary conditions this method 
allows one to recover the support function without a priori knowledge of the boundary 
condition. If the obstacle is not convex, the method recovers the convex hull of the 
obstacle. Numerically one can recover the obstacle, after its convex hull is found, by 
using Modified Rayleigh conjecture method, introduced in [TU], or by a parameter-fitting 
method. 

In ^S] there is a formula for finding an acoustically soft obstacle from the fixed- 
frequency scattering data. It is an open problem to develop an algorithm based on this 
formula. 

A numerical implementation of the Linear Sampling Method (LSM), suggested in (2], 
consists of solving a discretized version of 

[ A{a,p)gdp = (14) 

where A is the scattering amplitude at a fixed k > 0, is the unit circle, a E S^, and z 
is a point on M^. 

Let F = {Aai,Pj}, i = 1,...,N, j = 1,...,N be the square matrix formed by the 
measurements of the scattering amplitude for incoming, and outgoing directions, 
then the discretized version of (14) is 

Fg = f, (15) 



where the vector f is formed by 



see [T for details. 

Denote the Singular Value Decomposition of the far field operator hj F = USV^ . 
Let Sn be the singular values of F, p = U^i, and fi = V^i. Then the norm of the sought 
function g is given by 

^1 |2 

ii^ir = E^- (17) 

n=l " 

A different LSM is suggested by A. Kirsch in In it one solves 

(F*F)^/^g = f (18) 
instead of (15). The corresponding expression for the norm of Q is 

11^11^ = 5^^. (19) 

n=l 

A detailed numerical comparison of the two LSMs and the linearized tomographic inverse 
scattering is given in jl]. 

The conclusions of jT], as well as of our own numerical experiments are that the 
method of Kirsch (19) gives a better, but a comparable identification, than (17). The 
identification is significantly deteriorating if the scattering amplitude is available only 
for a limited aperture, or if the data are corrupted by noise. Also, the points with the 
sma//esi values of the ||^|| are the best in locating the inclusion, and not the largest one, 
as required by the theory in [S] and in P]. In Figures 1 and 2 the implementation of the 
Colton-Kirsch LSM (17) is denoted by gnck, and of the Kirsch method (19) by gnk. The 
Figures show a contour plot of the logarithm of the ||^||. The original obstacle consisted 
of two circles of radius 1.0 centered at the points {—d, 0.0) and {d, 0.0). The results 
of the identification for d = 2.0 are shown in Figure 1, and the results for d = 1.5 are 
shown in Figure 2. Note that the actual radius of the circles is 1.0, but it cannot be 
seen from the LSM identification. Also, one cannot determine the separation between 
the circles, nor their shapes. Still, the methods are fast, they locate the obstacles, 
and do not require any knowledge of the boundary conditions on the obstacle. The 
Support Function Method ([S], [7j) showed a better identification for the convex parts of 
obstacles. Its generalization for unknown boundary conditions is discussed in ^3]. The 
LSM identification was performed for the scattering amplitude of the obstacles computed 
by the Boundary Integral Equations method, see No noise was added to the synthetic 
data. In all the experiments we used k = 1.0, and N = 60. 
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gnk, k=1.0, d=2 . 




gnck, k=1.0, d=2 . 




Figure 1: Identification of two circles of radius 1.0 centered at {—d, 0.0) and {d, 0.0) for 
d = 2.0. 
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gnk, k=1.0, d=l . 5 




Figure 2: Identification of two circles of radius 1.0 centered at {—d, 0.0) and {d, 0.0) for 
d = 1.5. 
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